't Hooft-Polyakov Monopoles in an Antiferromagnetic Bose-Einstein Condensate 
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We show that an antiferromagnetic spin-1 Bose-Einstein condensate, which can for instance be 
created with ^^Na atoms in an optical trap, has not only singular line-like vortex excitations, but 
also allows for singular point-like topological excitations, i.e., 't Hooft-Polyakov monopoles. We 
discuss the static and dynamic properties of these monopoles. 
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Introduction. — Quantum magnetism plays an im- 
portant role in such diverse areas of physics as high- 
temperature superconductivity, quantum phase transi- 
tions and the quantum Hall effect. Moreover, it now 
appears that magnetic properties are also very impor- 
tant in another area, namely Bose-Einstein condensation 
in trapped atomic gases. This is due to two independent 
experimental developments. The first development is the 
realization of an optical trap for ^^Na atoms jl], whose 
operation no longer requires the gas to be doubly spin- 
polarized and has given rise to the creation of a spin-1 
Bose-Einstein antiferromagnet The second develop- 
ment is the creation of a two-component condensate of 
^^Rb atoms , which by means of rf-fields can be ma- 
nipulated so as to make the two components essentially 
equivalent Q|. As a result also a spin- 1/2 Bose-Einstein 
ferromagnet can now be studied in detail experimentally. 

The spin structure of these condensates has recently 
been worked out by a number of authors and also 
the first studies of the line-like vortex excitations have 
appeared An immediate question that comes 

to mind, however, is whether the spin degrees of free- 
dom lead also to other topological excitations that do not 
have an analogy in the case of a single component Bose- 
Einstein condensate. The answer to this question is in 
general affirmative. Indeed, we have recently shown that 
ferromagnetic Bose-Einstein condensates have long-lived 
skyrmion excitations, which are nonsingular but topo- 
logically nontrivial point-like spin textures |l^]. More- 
over, we show here that also spin-1 Bose-Einstein an- 
tiferromagnets have point-like topological excitations. 
In particular, there exist singular point-like spin tex- 
tures, which we call 't Hooft-Polyakov monopoles be- 
cause of their analogy with magnetic monopoles in par- 
ticle physics [|l^. Having done so, we then turn to the 
investigation of the precise texture and the dynamics of 
these monopoles. 

Topological considerations. — To find the topological 
excitations of a spin-1 Bose-Einstein condensate, we need 
to know the full symmetry of the macroscopic wave func- 
tion 'l'(r) = ■\/n(r)C(r), where n(r) is the total density of 
the gas, C(r) is a normalized spinor that determines the 
average local spin by means of (F)(r) = CUi")FC(r), and 



F are the usual spin matrices obeying the commutation 
relations [Fa^Ffj] = ieaPfF-^. Note that here, and in the 
following, summation over repeated indices is always im- 
plicitly implied. From the work of Ho § we know that 
in the antiferromagnetic case the mean-field interaction 
energy is minimized for (F)(r) = 0, which implies that 
the parameter space for the spinor C(r) is only x 5*^ 
because we are free to choose both its overall phase and 
the orientation of the spin quantization axis. Introducing 
the superfluid phase z9(r) and the unit vector field m(r), 
this topology can also be understood explicitly from the 
fact that all the spinors 



C(r) 



V2 



-TO:r(r) + imy{v) 
V2m^(r) 



= „»i?(r) 



CAF(r) (1) 



have a vanishing average spin and hence are locally de- 
generate. 

What does this tell us about the possible topolog- 
ical excitations? For line-like defects or vortices, we 
can assume ^(r) to be independent of one direction and 
the spinor represents a mapping from a two-dimensional 
plane into the order parameter space. If the vortex 
is singular this is visible on the boundary of the two- 
dimensional plane and we need to investigate the proper- 
ties of a continuous mapping from a circle into the or- 
der parameter space G, i.e., of the first homotopy group 
TTi{G) H. Since 7ri(S'i x S^) = Z, we conclude that 
an antiferromagnetic spin-1 condensate can have vor- 
tices with winding numbers that are an arbitrary integer. 
Physically, this means that by traversing the boundary 
of the plane, the spinor can wind around the order pa- 
rameter an arbitrary number of times. 

Similarly we can also discuss singular point-like de- 
fects. Since the boundary of a three-dimensional gas 
is the surface of a three-dimensional sphere, singular 
point-like defects are determined by the second homo- 
topy group 'K2{G) [|l3|. Because also 712(6'^ x S"^) = Z, 
such topological excitations thus indeed exist in the case 
of a spin-1 Bose gas with antiferromagnetic interactions. 
In view of the work of 't Hooft and Polyakov we refer to 
these excitations as monopoles, although it would also be 
justifiable to call them singular skyrmions. In contrast 
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to the nonsingular skyrmions in the Bose-Einstein fer- 
romagnets, which inherently are nonequilibrium objects, 
the 't Hooft-Polyakov monopoles turn out to be thermo- 
dynamically stable excitations as we show next. 

Monopole texture. — The grand-canonical energy of 
the spinor condensate can be obtained from the usual 
Gross-Pitaevskii theory, which for the restricted param- 
eter space given in Eq. (|^), leads to the expression 

E[n, C]= Jdr V/(r) (-^^ + ^trap(r) - ^l 

+ ylV'(r)P + ^[Vm(r)]^)^(r) , (2) 

where i^{r) — n(r)e''''^''^ is the superfluid order param- 
eter, T4rap(r) = muP'Y^ /2 is an isotropic harmonic trap- 
ping potential, /x is the chemical potential, and (/„ = 
iiranh^/m is the appropriate coupling constant for den- 
sity fluctuations. Minimization of this energy determines 
both the spin texture m(r) and the density profile n(r) 
of the 't Hooft-Polyakov monopole jlj] . Since gradients 
in the spin texture do not couple to the superfluid phase 
i?(r), we can at this point already conclude that the pres- 
ence of a monopole will not induce any superfluid flow in 
the atomic cloud. From now on we therefore no longer 
consider this degree of freedom. 

Interestingly, the spin texture is uniquely determined 
by the fact that it should have a topological winding num- 
ber (l| 

W ^ [ dr eijk€a0jdimadjmpdkm^ (3) 
Svr J 

equal to 1 and that it should also minimize the gradient 
energy 

i?grad[n,C] -Jdr !^[Vm(r)]2 . (4) 

As we discuss in more detail in a moment, the latter re- 
quires the spin texture to be as symmetric as possible. In 
combination with the first requirement, we thus conclude 
that mHp(r) = r/r for a monopole in the center of the 
trap. Indeed, a spin texture with the same winding num- 
ber can be obtained by rotating mHp(r) by an arbitrary 
rotation matrix R(r) that only depends on the radial dis- 
tance r. As a result the gradient energy is turned into 
the suggestive form 

i?grad[n,C] = [dr !^|[V -zA(r)]mHp(r)p , (5) 

which brings out even more clearly the analogy with the 
0(3) gauge theory studied by 't Hooft and Polyakov. 
In our case the vector potential is, however, defined by 
A(r) = i(R^-'(r)VR(r)) and is always orthogonal to 
the gradient of mnp (r) . The gradient energy is therefore 
minimized for R(r) = 1, as anticipated previously. 
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FIG. 1. Spin texture and density profile of the 't 
Hooft-Polyakov monopole. The solid line shows the exact 
numerical result, whereas the dashed line shows the analytic 
long distance behaviour n(r) ~ {fi/gn){l — 2(^/r)'^). 

Substituting our hedgehog solution for the spin tex- 
ture into Eq. (|^), wc see that the density profile of the 
spinor condensate is found from a Gross-Pitaevskii equa- 
tion with a centrifugal barrier equal to fi^ /nir"^. Consid- 
ering first the homogeneous case and writing the super- 
fluid order parameter as ip{r) — [fJ-fgny^^ f{r/^), with 
^ — {h^ /2mfi)^^^ the correlation length, we obtain ex- 
plicitly that 

This equation cannot be solved analytically, but its nu- 
merical solution is shown in Fig. 1. What is most 
important for our purposes, however, is the large dis- 
tance behavior of the density profile. Neglecting the 
gradient terms in the left-hand side, we easily find that 
f{p) ~ 1 — for /9 ^ 1. Furthermore, the monopole 
is clearly seen to possess a core with a typical size of 
the order of the correlation length. Inside this core the 
density is strongly reduced, which offers the opportunity 
to detect monopoles by the same expansion experiments 
that have recently also been used to observe vortices . 

Having obtained the spin texture and the density pro- 
file of the 't Hooft-Polyakov monopole, we are now in 
a position to determine also its energy. Placing the 
monopole in the center of a spherical volume with a large 
radius i? S> ^, we obtain 

^^HP — —{R~ i?corc) , (7) 

where i?corc is the effective core size of the monopole. Its 
calculation requires knowledge of the complete density 
profile and we find numerically that i?coic — 1.4^. The 
monopole energy thus diverges linearly with the system 
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size, which imphes that in the thermodynamic hmit only 
pairs of monopoles with opposite winding numbers re- 
quire a finite energy for their creation. 

The physical interpretation of Eq. (|^) is that to calcu- 
late the monopole energy we only need to evaluate the 
gradient energy for a fixed density equal to /i/(?„, but 
restrict the integration to the volume outside a spheri- 
cal core region with radius i?core- This interpretation is 
particularly useful for a trapped spinor condensate in the 
Thomas-Fermi limit, when the size i?TF = 
of the condensate is much larger than the correlation 
length ^. For a monopole at position u near the cen- 
ter of the trap, we find in this manner that 



Shp[u] 



2fj. 
3a„ 



i?TF 1 



3u^ 



3u4 



(8) 



if we neglect the core contributions that are smaller by a 
factor of C/i?TF- The significance of this result will be- 
come clear once we understand the dynamical properties 
of the 't Hooft-Polyakov monopole. 

Monopole dynamics. — In first instance we expect the 
dynamics of the monopole to be determined by the action 
^[n, C] = / dt{T[n,(] ~ E\n.,Q\) with a time-derivative 
term that is equal to 



r[n,C]= y" dv n(r,t)C^(r,t)i?i^C(r,i) 



(9) 



in the Gross-Pitaevskii theory. However, when we re- 
strict ourselves to the antiferromagnetic spinor C^af (r, t) , 
the time-derivative term in the action exactly vanishes 
due to the normalization condition m^(r, = 1. To find 
any dynamics for the monopole we thus need to consider 
also fluctuations that bring the spinor condensate out 
of the antiferromagnetic order parameter space. In the 
Thomas-Fermi limit, the relevant dynamical part of the 
action thus becomes 

= 1 dt(T[nX]-J dr |[n(r,i)(F)(r,i)]2) , (10) 

with Qs = Anash^ /m the appropriate coupling constant 
for spin-density fluctuations. Using this action we can 
now investigate the effect of the above mentioned fluc- 
tuations by substituting C(r,t) = CAF(r,i) +SC{r,t) and 
expanding the action up to quadratic order in 5C,{T,t). 
Solving then the Euler-Lagrange equation for S({r, t) and 
substituting the solution back into the action, we ulti- 
mately find the desired low-frequency result [[l7[ 



'5'dyn[", (]= J dt j dr ^ 



9m(r, t) 
Ft 



(11) 



While performing the calculation, we must make sure 
that we are not considering fluctuations of the spinor 



within the antiferromagnetic order parameter space. 
This requires the matrix elements (l^-p{r,t)F6({r,t) to 
be real, because for fluctuations within the antiferro- 
magnetic parameter space we have in lowest order that 
= (F)(r,t) = ClF(r,t)F,5C(r,i) +,5Cnr,i)FCAF(r,t). 
Moreover, the normalization of the spinor requires also 
that ClF(i',i)'5C(r,t) =0. 

The importance of this result is twofold. First, from 
Eqs. (|^) and ( pT| ) we sec that at the quantum level the 
dynamics of the space-independent part of the vector 
field m(r) is governed by the following time-dependent 
Schrodinger equation 



(12) 



for the wave function \E'(m, i). It thus corresponds ex- 
actly to a quantum rotor with a moment of inertia equal 
to / = h^Voip)/ gs, where Vb(/i) is the total volume of the 
spinor condensate in the Thomas- Fermi limit. In an har- 
monic trap the moment of inertia is thus proportional 
to the 3/5 power of the total number of atoms. The 
eigenstates of this Schrodinger equation are the spherical 
harmonics Ygjv/3(m). In this way we thus recover the 
fact that according to quantum mechanics both the to- 
tal spin of the Bose-Einstein antiferromagnet as well as 
its projection on the quantization axis must always be 
an integer. More precisely, since the ground state wave 
function is given by Fo.olm) = l/\/47r, we have actually 
shown that the many-body wave function of the antifer- 
romagnetic spinor condensate is a singlet state exactly 
Note that physically this phenomenon is equiva- 
lent to the way in which 'diffusion' of the overall phase of 
a Bose-Einstein condensate leads to the conservation of 
particle number (l^Jl^ . The main difference is that here 
the 'diffusion' takes place on the surface of a unit sphere 
instead of on a unit circle. 

Second, and most important for our purposes, we can 
now determine the single monopole dynamics, by using 
the ansatz m(r,t) — mHp(r — u(t)) for the texture of a 
moving monopole, which is expected to be accurate for 
small velocities du{t)/dt and, in the inhomogeneous case, 
near the center of the trap where u/R^p <C 1. Substi- 
tuting this ansatz into Eq. (11) and remembering also 
Eq. (^), we find that the action for the center-of-mass 
motion of the monopole becomes precisely that of a mas- 
sive particle 



(13) 



with a mass given by tohp — 2mi?TF/as- The semi- 
classical equation of motion for the monopole position is 
therefore simply Newton's equation 



mnp- 



d^ujt) 
dt^ 



(14) 
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In view of the fact that the energy of the 't Hooft- 
Polyakov monopole decreases as the distance to the cen- 
ter of the trap increases, we conclude that in general the 
monopole is always accelerated to the boundary of the 
condensate. Typically it will reach that boundary in a 
time interval of order 



' mR^p an 



2^1 




(15) 



which for present-day experiments is sufficiently long to 
be able to observe the monopole once it is created Q . In 
this context it should be noted that all our calculations 
are performed at zero temperature. In the presence of a 
normal component, the monopole experiences damping, 
which increases the time needed to reach the edge of the 
spinor condensate. 

Summarizing, we have investigated the most important 
equilibrium and nonequilibrium properties of a 't Hooft- 
Polyakov monopole in a trapped, antiferromagnetic Bose- 
Einstein condensate. A further direction of research is 
an ah initio calculation of the above mentioned friction 
force on the monopole. In fact, we expect the thermal 
cloud to lead not only to dissipation but also to noise. 
Both effects can be conveniently treated within the gen- 
eral framework of the stochastic field theory that was 
developed previously for nonequilibrium phenomena in 
partially Bose-Einstein condensed gases . Another in- 
teresting topic is the interaction between two monopoles 
and the many-body properties of a gas of these topologi- 
cal objects. Finally, a very important experimental prob- 
lem is the creation of a 't Hooft-Polyakov monopole. Of 
course, pairs of monopoles with opposite winding num- 
bers can in principle be created in a thermal quench or by 
sufficiently shaking up the spinor condensate. However, 
a more controled creation mechanism is desirable. There- 
fore, we are presently exploring if a single monopole can 
be created by an appropriate tailoring of the detuning, 
the intensity and the polarization of two pulsed Raman 
lasers. 
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